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SUMMARY 


A study of some of the important aerodynamic factors affecting the 
directional stability of supersonic airplanes is presented. The mutual 
interference fields between the body, the lifting surfaces, and the sta- 
surfaces are analyzed in detail. Evaluation of these interfer- 
ence fields on an approximate theoretical basis leads to a method for 
predicting directional stability of supersonic airplanes. Body shape, 
wing position and plan form, vertical tail position and plan form, and 
ventral fins are taken into account. Estimates of the effects of these 
factors are in fair agreement with experiment. 


INTRODUCTION 


Aircraft designed for supersonic flight have problems of directional 
instability. The difficulty arises primarily because the center of grav- 
ity of the supersonic fuselage -engine combination is well aft of the 
fuselage nose. The resulting "short" afterbody then restricts the "moment 
arm available to the stabilizing surfaces. Current methods of compensat- 
ing for this moment am penalty are largely limited to the use of large 
vertical tails and auxiliary fins, such as "strakes" or ventral fins. 
Destabilizing aerodynamic effects which compound the difficulty of main- 
taining directional stability occur with increasing angle of attack and 
Mach number. 


The important aerodynamic factors contributing to directional insta- 
bility of supersonic airplanes were briefly pointed out in reference 1. 
These factors are reviewed here to outline the nature of the study which 
will subsequently be undertaken in this report. The directional charac- 
os of an airplane are dominated by those of the body and the ver- 
tical tail. Insight into the directional stability characteristics of a 
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supersonic airplane is then provided by considering the properties of a 
body-tail configuration. Vertical -tail surfaces are used to compensate 
for the directional instability of the fuselage and to give positive 
directional stability to the combination. The stabilizing moment of the 
combination about the fuselage center of gravity can be separated into 
moments caused by the tail and the body alone. The unstable moment coef- 
ficient of the fuselage is nearly independent of Mach number. The stabi- 
lizing moment contributed by the tail decreases with increasing Mach 

number approximately as and the directional moment coefficient 

of the combination then also decreases with Mach numbei , as shown in 
sketch (al). Further decreases in the moment coefficient occur when the 
airplane is at combined sideslip and large angle of attack. This addi- 
tional loss in stabilizing moment is caused by destabilizing elements in 
the fuselage sidewash wake acting on the vertical -tail panels. This wake 
is generally characterized by the presence of two body vortices whose 
strength increases approximately in proportion to the square of angle of 
attack. At angles of sideslip, the vortex pair is asymmetrically disposed 
with respect to the tail plane and causes a destabilizing sidewash with, 
the attendant stability characteristics shown in sketch (a2). The combined 



Sketch (a) 


effect of Mach number and angle of attack in many instances can cause an 
airplane to become directionally unstable. 
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It is the purpose of this report to present a method for quantita- 
tively predicting the foregoing destabilizing effects along with the less 
obvious aerodynamic effects of a wing on directional stability of a body 
and body-tail in combination. 


SYMBOLS 


The results are presented as standard NACA coefficients of forces 
and moments. The data are referred to the body axis system with the 
reference center of moments located as indicated in figure 1. 
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body vertical semiaxis , in. 
body horizontal semiaxis , in. 
body 


crossflow drag coefficient 


wing root chord at body juncture , in. 


distance to center of pressure from the juncture of wing leading 
edge and body, in. 


yawing -moment coef f i ci ent , 


yawing moment 


c 3 00 W 2s W 
ide force 




side-force coefficient, 
body diameter, in. 
horizontal stabilizer 
vortex interference coefficient 

body-w T ing interference coefficient due to angle of attack 
interference coefficient 
body length, in. 

distance of center of pressure from body nose, in. 

distance to moment reference point from body nose, in. 

length of body ahead of juncture of the wing leading edge and 
body, in. 
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M local Mach number 

M c crossflow Mach number, M^sin \J a^+p 2 

free-stream Mach number 
N yawing moment, lb in. 

p local static pressure, lb/sq in. 

p^ free-stream static pressure, lb/sq In. 

q local dynamic pressure, lb/sq in. 

free-stream dynamic pressure, lb/sq in. 
r body radius, in. 

Sjj ventral-fin span from body axis, in. 

sy vertical -tail span from body axis, in. 

wing semispan, in. 

S-g maximum cross-sectional area of body, sq in. 

Sp projected area of body on vertical plane, sq in. 

Sty total area of wing including portion blanketed by body, sq in. 

U ventral fin 

v lateral component of free-stream velocity 

V upper vertical tail 

w vertical component of free-stream velocity 

W wing 

x body coordinate measured from nose 

X streamwise coordinate 

Z vertical coordinate 

a angle of attack, radians unless noted 


3 


angle of sideslip, radians unless noted 
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ratio of specific heats , 1*40 
Body vortex strength 
index denoting increment 

angle Between tail leading edge and Body center line, deg 

corrective coefficient to Body-wing interference coefficient of 
slender-Body theory 

effectiveness factor involved in wing expansion or compression 
fields 

proportionality constant 
taper ratio 


Subscripts 


Body 

horizontal tail 
vortex separation 
ventral fin 
upper vertical tail 
wing 

combination of Body and upper vertical tail 
combination of Body, upper vertical tail, and ventral fin 
combination of Body and wing 

combination of Body, wing, and upper vertical tail 

combination of Body, wing, upper vertical tail, and horizontal 
tail 

combination of Body, wing, upper vertical tail, and ventral fin 

in presence of (e.g., B(V), Body in presence of upper vertical 
tail) 
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3 coefficient derivative (e.g., ( )p, derivative of coefficient in 

parentheses with respect to sideslip angle p) 

cp cross -coupling coefficient 


Superscripts 


4 - compression field 

expansion field 
F sweptforward wing 

R reversed wing 


ANALYSIS 


As a guide to development of the estimation procedures of the 
present report , a systematic series of tests on a variety of airplanes 
was required. In particular, tests were nee led to reveal the directional 
stability increments contributed by each component in a progressive 
build-up from a basic body to a complete airplane. The bulk of such 
experimental data were provided through tests of the configurations shown 
in figure 1 in the Ames 1- by 3 -foot supersonic wind tunnel No. 2 at Mach 
numbers 1.97 to 2.94. A description of this facility is given in refer- 
ence 2. Additional data were taken from references 3 to 5. 

The estimation procedures to be presented in the present report are 
an extension to those of reference 6. The linear theory is reviewed for 
completeness of presentation, but emphasis Is given to nonlinear effects 
attending large angles of inclination. 

To determine the aerodynamic characteristics of a body-component 
combination the body and component characteristics must be known inde- 
pendently, as well as the mutual interactions existing between the body 
and components in combination. The method of determining mutual inter- 
ference effects at low angles is based on the result from slender-body 
theory that the ratio of the change in side force due to the addition of 
a surface to a body to the side force developed by the surface alone is 
in the same ratio as the change in the apparent mass (of the body-surface 
cross section in sideslip) due to the addition of the surface to the 
apparent mass of the surface alone. The change or increment in apparent 
mass due to the addition of a component to a pre-existing configuration 
causes a redistribution of force interactions among the components. The 
distribution of this incremental side force among the components is not 
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determined in the analysis to he presented; only the total change is 
known. This incremental side force is simply charged to the component 
that produces the change. As a consequence of this procedure , the order 
in which a configuration is constructed affects the ''interference coef- 
ficients." Accordingly, the interference coefficients presented must he 
used in the order prescribed. Two sets of interference coefficients are 
presented: (l) coefficients which are applicable to uncoupled wing -tail 

configurations (aircraft having no interaction between the wing and 
directional stabilizers), and (2) coefficients applying to closely coupled 
wing-tail configurations (aircraft with wing in close proximity to the 
directional stabilizers). Symbolic designations of the incremental inter- 
ference coefficients K or apparent mass ratio changes along with associ- 
ated force designations are illustrated in figure 2 for coupled and 
uncoupled wing-tail configurations. The symbols are arranged in a sequence 
indicating the order of configuration build-up to be employed. Isolated 
component forces, denoted Y, are determined by estimate or experiment. 
Figures giving values for the interference coefficients K are introduced 
at appropriate sections in the text. The foregoing is essentially a review 
of the method of reference 6 for estimating the directional stability 
coefficients of airplanes at low angles of sideslip and attack. 

In the following sections the theoretical directional stability 
coefficients of various body- component combinations are derived, beginning 
with an initial consideration of the directional characteristics of four 
different bodies. Methods applicable to determining the incremental 
directional coefficients resulting from the addition of various components 
are then considered at each stage of a configuration build-up. Linear 
theory applicable at low angles is presented for each case, followed by 
methods for estimating the nonlinear effects introduced by large angles 
of inclination. Correlation of experiment and theory at both low and 
high angle ranges is discussed. 


Body Alone 


Pure sideslip.- Available theories generally give satisfactory values 
for aerodynamic forces on a body. On the other hand, simple methods of 
estimating the center of pressure of these forces are unavailable. In 
view of the lack of accurate theory for determining body centers of pres- 
sure, estimates of body-alone yawning moments are omitted in the present 
report. Consideration, however, Is given to determining side force of 
bodies for the sake of maximum completeness of presentation. The value 
of the side -force coefficient near (3=0° for circular bodies is given 
by slender-body theory as 
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L B 


20 


(la) 
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based on maximum body cross-sectional area. At large angles of sideslip, 
force coefficients of bodies become nonlinear with angle of inclination 
as a result of the effects of viscous separation and are then grossly 
underestimated by slender-body theory. The nonlinear increase in the 
force coefficient can be empirically taken into account, by the method 
of reference 7, through an additive body crossflow coefficient term 



where cq c is the two-dimensional drag coefficient for the body cross 

section and Sp is the projected body area normal to the crossflow direc- 
tion. It is assumed that, at supersonic speeds, the crossflow about the 
body at a given longitudinal location is unaffected by "end effects" of 
body nose or base. The use of two-dimensional crossflow coefficients is 
then possible. An investigation of bodies having elliptical cross sec- 
tions (refs. 2 and 8) shows that the experimental side -force coefficient 
slopes are predicted with fair accuracy wnen equation (la) is modified 
by the ratio of the semiaxes lengths, a/b, in accordance with slender- 
body theory. The investigation shows further that the force coefficients 
of an elliptical and circular body of the same base reference area are 
in the ratio a/b over a wide range of sideslip angles. Thus, the term 
a/b affects the potential (linear) and viscous (nonlinear) terms of body 
force coefficients in the same degree. In most cases bodies are either 
nearly circular or elliptical In cross section. An approximate method 
of determining the side-force coefficient for these cases is based on two 
assumptions: (l) The linear term is given by slender-body theory and is 



(2a) 


where a/b pertains to the semiaxis ratio at the maximum cross-sectional 
area, Sp, location. (2) Essentially two-dimensional crossflow prevails 
at each section and the local crossflow coefficient is given, by general- 
izing the results of reference 8, as the local semiaxes ratio a/b as a 
function of x times the crossflow coefficient of a circular cross sec- 
tion. The nonlinear crossflow force is then the integrated effect of 
local crossflow coefficients over the projected body area. Within the 
limitation of these assumptions the total side-force coefficient is 
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Equation (2b) gives side-force coefficierts that are in fair accord with 
experimental values, as will be shown. 
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Combined angles of attack and sideslip .- An expression giving the 
side-force coefficients for a circular cross-section body (a/b = l) in 
combined sideslip and angle of attack is given simply as the sideslip (p) 
component of the fo rce co efficient associated with the total angle of 
inclination (a 1 = J a 2 +p 2 ) ; that is, 


S3 

C Y = ~ 2a' 



+ 


C *c 




= ^ 2p + c dc J p/a 2 + p 2 (3) 

For "bodies of noncircular cross section (a/b / l), “the force vector does 
not necessarily lie in the plane of maximum body inclination a 1 . This 
fact precludes simple estimates of side -force coefficients of other than 
circular cross-section bodies at combined pitch and sideslip angles . 1 

Comparison of estimate and experiment . - In order to evaluate the 
side -force coefficients of bodies, the crossflow drag coefficient for 
circular bodies must be known. At low, subsonic crossflow velocities, 
the drag coefficient is 1.2. At high velocities (M c > 0.4), c& c 

varies significantly and the variation must be taken into account. Exper- 
imental values for at the Reynolds number of the present test were 

abstracted from reference 10 and are plotted as a function of crossflow 
Mach number, M c , in figure 3. These estimates for c^ c were used in 
evaluating equations (2b) and (3). Figure 4 is a plot of experimental and 
estimated values of Cy as a function of sideslip angle 3 for three 
body models of the present test at = 2.94. Good correlation between 
experiment and estimate is shown over the angle range investigated. It 
can be concluded that, for the range of body shapes investigated, equa- 
tion ( 2 b) adequately predicts the effect of "shape" on body side -force 
coefficients . 

In figure 5 experimental and estimated values for the side-force 
coefficients of circular body models are plotted as a function of angle 
of attack for two sideslip angles . The estimated values are given by 
equation ( 3 ). The effect of angle of attack on the side-force coefficient 
is fairly well predicted. A point of particular interest lies in the 
fact that the effect of increased body length on the side -force coeffi- 
cient (the difference of Cy between bodies B x and at (3 = 5 °) is 
correctly predicted. This prediction of side -force change due to body- 
length change has an important bearing on the determination of the effects 
of body-wing interaction and will be discussed in a later section. 

■‘"Consideration of other than circular cross-section bodies is given 
in reference 9 . 
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Body-Tail Configuration 


In the present method, it is assumed that the yawing -moment coeffi- 
cients of the body alone have been experimentally determined. These 
coefficients are used in calculating values for the directional coeffi- 
cients of the body in combination with wings and stabilizing surfaces. 

Body-tail configuration in pure sideslip .- The total side force of 
the combination of a body with a tail is the sum oi the forces contributed 
by the body alone and by the tail in presence of the body. This sum is 
expressed in figure 2 for a body, B, and upper vertical tail, V, as 

Y BV = y B +k V(B) y V 

It is convenient, in comparisons with experiment, to express equations 
for the directional coefficients in a form which denotes the stabilizing 
contribution of the added component. Accordingly, equation (4a), expressed 
in coefficient form as the contribution of the tail In the presence of 
the body, is written 

Cy bv' Cy b = k v(b) c y v 


The side-force coefficient Cy v is that of ;he upper tail panel when 
mounted on an infinite reflection plane (equivalent to one-half the coef- 
ficient of two such panels joined together). The side-force coefficient 
Cy B is that of the body alone. Values of C V ( B ) abstracted from refer- 
ence 6 are indicated by the solid curve in figure 6 as a function of ratio 
of body radius to tail span. 

The K coefficients take into account nutual interaction between 
components; for example, Ky( B ) includes the effect of interference on V 
due to B and the effect of force carryover onto B due to V. This is 
expressed symbolically as ^j(B) ~ ^■V+B- ) -V + ^-V tB • ^ breakdown of tnese 

interference effects is not essential In determining total side force; 
however, the division of loading between body and tail and the correspond- 
ing centers of pressure (T) must be known in order to determine yawing 
moment accurately. Since the present method does not give a division of 
the interference forces, an approximate method of determining moment 
values Is employed. This method is developei as follows: 

N BV = N B + %+B-V +W V -B 
or 

it \ / 4y+B->V _i M v , Z V-KB- Z M v 'N 

%v = Nb+Ov-^m) ( — = — - — y v+b-»v + -= — - — y V-*B ) 

\ ly-i M l Y ~ l M / 


(5a) 
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The distance to the tail center of pressure from the moment reference 
point (iy-Zjyj) in the absence of body interference is introduced in equa- 
tion (5a). This value is not significantly different from that of the 
tail in the presence of the body ( W+B-»V"^m) • is assumed that the 

center-of -pressure location of the force distribution on the body due to the 
tail is approximately that of the isolated tail. Within these 

approximations, and with Yy( B ) = Yy +B _*y4Yy_> B = Ky( B )Yy, equation (5a) 
can be simplified and written in coefficient form as the yawing -moment 
contribution of a tail in the presence of a body 

Cn BV -Cn B = Ky( B )C Yv (5b) 

The side-force-coefficient increment due to a tail assembly including 
a ventral fin is 

C Y bvu - C Y b = K V(B) c Yy+ K U(BV) c Y u (6) 

Values for Ky^y) are also given in figure 6 by the dashed line param- 
eters which take into account the effect of the upper tail. The definition 
of Cy^ is similar to that of Cy^. The incremental yawing-moment- 
coefficient equation is developed through steps similar to those leading 
to equation (5b) to be 

c n B yu- C n B = (i 5S~) K V(B)C Y y + (- ) K U(BV) C Yu (7) 

Body-tail in combined sideslip and angle of attack . - The effects of 
angle of attack on the directional coefficients of body-tail configura- 
tions , and of configurations in general, arise from two distinct aerody- 
namic interference phenomena. These are (l) a cross -coupling effect of 
sideslip and angle of attack induced crossflow velocities, discussed in 
reference 11, and (2) an induced si dewash due to vortices (particularly 
those forming about the body at high angles of attack). These two effects 
are discussed in turn. 

The cross -coupling interference phenomenon is taken into account by 
the factor K^, developed in reference 11, with which the modifying effect 
of angle of attack on the tail contribution to the side-force coefficient 
of a body-tail combination can be determined. The increment in side-force 
coefficient is given by the equation 

%(3,a) - % 


( 8 ) 
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The factor is plotted as a function of ratio of panel semispan to 

body radius, r/sy, in figure 7 for a planar and a cruciform body-tail 
configuration. It is to be noted that values are given only for equispan 
double and cruciform panel body configurations. Solutions for a single 
panel and body combination are not available at present. However, the 
experimental results of reference 11 show that panel -to -panel effects are 
negligible for a planar configuration with r/sy = 0.273 and M^ = 2. It 
is then indicated that for Mach numbers of 2 or greater, the K<p applying 
to planar (two surface) configurations can also be applied to a single 
surface configuration. In most cases considered in this report the Mach 
number is 2 or greater and the values presented in figure 7 are adopted. 

In order to estimate the effects of body-vortex induced sidewash on 
the tail side-force coefficient, the position of the vortices with respect 
to the tail plane and the strength of the vortices must be known. Theo- 
retical methods for determining the position and strength of body vortices 
involve very laborious calculations even for circular bodies and are 
unavailable for noncircular bodies. An experimental correlation of these 
quantities is presented in reference 12 for circular bodies and provides 
a simple estimate for vortex strengths and positions in lieu of theory. 
With these quantities in hand, the effect of body vortices on the upper 
tail side force is related t o the t ail -alone side force developed at the 
angle of inclination a' = si a 2 +0 2 by the factor Aiy and the body vortex 
strength f B 


CY v(r B ) 


= Ai-v 


r B 


2n si v 2 +w 2 r 



CY V(a' ) 


(9) 


Values for Aiy are determined, as shown by the example in the appendix, 
from figure 8. Interference factors (iy) for this chart were determined 
by the equations of reference 13 for the effect of a single external vor- 
tex and its body image on an upper tail pan si. (The effect of two exter- 
nal vortices is found by superposition. ) rhis chart applies strictly to 
a tail panel of taper ratio A = l/2 and r/sy = 0.2 corresponding approxi- 
mately to that of model BjjVs. The use of this chart for panels of other 
taper ratios will be discussed in a following section on experimental 
considerations . 

The side-force coefficient contributicn of an upper tail to a body- 
tail combination at combined sideslip and angle of attack is then given 
by the sum of equations (4), (8), and (9) 

= ( Kv ( b ) " ban Cyy CY V +CY V(r B ) 


c Y B v“ C Y B 


( 10 ) 
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The yawing -moment coefficient is 


C n B v‘% 



Kcpa \ 

tan CY V +CY v(r B ) 


( 11 ) 


The directional coefficients for a body-tail configuration including a 
ventral fin are obtained in a similar manner but with the cross -coupling 
term additive. 


and 


c y bvu - Cy b 


k v(b) - 


Kcpa 
tan e 


V- 


C Yv + 


k U(BV) 


Kcpa 
tan € 


UJ 


° Y U + 


: *v(r B ) + %(i V ) 


' B ‘ 


( 12 ) 


c n “ c n 
n BVU n B 


ly-l m 


2 s 


W 


2s w 


k v(b) - 


K U(BV) 


ify 


a 


tan e- 


K<p 


■V, 


a 


Cy +Cy . . 

Y v Y v(r B ) 


tan ejjJ Cy u +Cy U(i’ b ) 
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Comparison of estimate and experiment . - Estimated and experimental 
directional coefficients of body-tail combinations are compared on the 
basis of the coefficient values of the body and tail in combination minus 
that of the body alone. These comparisons are presented in figure 9 * 
Discussion is expedited by being directed primarily to the correlation 
between experimental and estimated side-force coefficients because the 
degree of correlation of yawing-moment data with estimates closely 
parallels that of the side force. 

The estimated values given apply strictly to the circular body-tail 
model B1V5 (fig. 9(a)). It is evident from the data that a strong desta 
bilizing effect occurs with increasing angle of attack. The estimated 
side-force coefficient, given by equation (4), is applicable at a = 0° 
(pure sideslip) and is in good accord with experimental values. If cross 
coupling is taken into account in the method (eqs. (4) and (8)), an appre 
ciable decrease in the directional coefficients results which, however, 
accounts only partially for the experimental decrease. Finally, equa- 
tion (10), which takes into account both cross coupling and body-vortex 
effects, is in good accord with experimental values over the angle-of- 
attack range of the test. Most of the remaining correlation figures, 

9 (b) to (f), show good agreement of estimate with experiment. It should 
be noted that figures 9( e ) and (f) give the incremental coefficients in 
derivative form at (3 ~ 0°. 
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The effect of body shape on the tail contribution to directional 
coefficients of body-tail combinations cannot be predicted quantitatively 
since the theoretical effects of body shape on vorticity and cross coupling 
are unknown. Accordingly, the effects of body shape were evaluated exper- 
imentally. To this end, coefficients for all three basic models are 
plotted in figure 9(a) to facilitate comparison. At the sideslip angle 
P = 5.1°, the coefficients of the tail in presence of both the elliptic 
body B 2 and the double elliptic body B 3 are less affected by angle of 
attack than in the presence of the circular body Bi. At sideslip angle 
3 = 10°, there is less effect of body shape among these models. 

A comparison of the side-force contribution of a tail on the extended 
nose circular model BjX with that of the same tail (V 2 ) on a rectangular 
afterbody model Ba (fig. 9( c )) indicates large influence of body shape. 

The rectangular model is considerably less stable over the angle -of -attack 
range than is the extended-nose model. A direct comparison of the effect 
of the two bodies on directional coefficients is not strictly justified 
since Ba has somewhat greater volume than and the strength of its 

vortex wake may be expected to be the larger. More important with regard 
to tail "effectiveness” are the different distributions of vorticity 
brought about by the difference in body shape of the two bodies. Vapor- 
screen pictures revealing this vorticity about bodies Ba and BiL in the 
vicinity of the tail surface are shown in f:.gure 10. Body Ba has three 
distinct regions of vorticity, the strongest lying in close proximity to 
the plane of the vertical tail. Body has two equal -strength vor- 

tices, one of which again lies in close proximity to the tail plane. The 
degree of destabilizing effect of body vertices on the vertical -tail sur- 
face is dependent on the vortex positions with respect to one another and 
on their relative position with respect to she tail plane. Vortices of 
equal strength (with opposite rotation) in :lose proximity to each other 
relative to their mean distance from the ta:.l surface induce small desta- 
bilizing sidewash. Vortices in close proximity to the body have less 
effect than when outboard. (The foregoing statements are readily verified 
by inspection of iy values in fig. 8.) Examination of the vapor-screen 
photographs in figure 10 show that in the case of body Ba the vortices 
are more widely separated from one another nhan are those about body BiL. 
Furthermore, the dominant (and destabilizing) vortex at the plane of the 
tail is at a greater tail spanwise distance from the body Ba than is 
the destabilizing vortex of body BiL- These differences in position are 
the probable cause for the relatively poor stability characteristics of 
BaV 2 in comparison with those of BipV 2 . It is to be noted that the 
cross-sectional shape of the afterbody of Ba (fig* 1(a)) has relatively 
sharp corners. A plausible illustration of how these "corners" can pro- 
mote the formation of vorticity into the complex wake is shown in fig- 
ure 10. In effect, corner vortices generated at sideslip and at angle of 
attack are swept in the direction of resultant flow to form the asymmetri- 
cal pattern shown. In view of this, deep angular body cross-sectional 
shapes may be undesirable since they can augment the unfavorable disposi- 
tion of body vorticity about the tail. 
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In determining the effects of vortex-induced sidewash on tail 
effectiveness, no account was taken of tail plan form as a variable. As 
stated previously, the influence of body vortices on a tail plan form 
corresponding to that of V5 was used as a basis for the chart of fig- 
ure 8- The general, though approximate, applicability of the chart for 
all plan forms (with r/sy - 0.2) is justified by experimental considera- 
tions. As in reference 13, it is assumed that the loss in tail loading 
at a given span location is proportional to the product of local chord 
length and the local change in flow inclination due to vortex-induced 
crossflow. The total tail -load loss is then determined by "strip" inte- 
gration of this loading loss over the exposed tail span. It is assumed 
in this procedure that a constant flow inclination with respect to the 
tail surface produces a span load distribution proportional to chord 
length distribution. In actuality, plan forms at moderate Mach numbers 
and low angles of inclination tend to develop a loading intermediate 
between elliptical and plan formlike in s panwise distribution. This tend- 
ency Is illustrated in sketch (b) in which comparison is made between 



Sketch (b) 

elliptical loading and normalized experimental loading for surfaces of 
rectangular and triangular plan form. It is evident that the assumption 
of plan formlike distribution of loading for each surface is not essen- 
tially any more valid than assigning an elliptical loading to both sur- 
faces. Since the actual loading does not correspond to either type, a 
straight-line (trapezoidal) loading was adopted as a reasonable compromise 
for the sake of computational simplicity. This trapezoidal loading was 
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used to determine the values presented in the chart of figure 8- This 
loading corresponds to the span distribution of chord length for the V5 
tail plan form. 

The foregoing concession made to computational simplicity is experi- 
mentally justified since the degree of correlation between experiment and 
theory does not vary significantly among the trapezoidal (fig. 9(a) ), tri- 
angular (fig. 9(t>))j and near rectangular (fig. 9(d)) plan-form tails 
investigated. In the case of the clipped tail (fig. 9 (d)), additional 
computations using the equations of reference 13 were required since body 
radius to span ratio was considerably larger than 0.2, precluding the use 
of figure 8. It is evident in figure 9 (d) that the ratio of stabilizer 
effectiveness at high angle of attack to its value of a = 0° is substan- 
tially reduced for the short or clipped span tail V 6 q in comparison to 
that of the unclipped tail V 6 . This is because a short span tail has a 
large portion of its surface exposed to the region of adverse sidewash 
existing between the destabilizing vortex core and the body surface. 

Additional correlation of estimated values with experiment is 
presented in figures 9( e ) and 9(f) for directional coefficient derivatives 
with (3 - 0° of two configurations, one of wnich includes a ventral fin, 
at Mach numbers of l.kl and 2 . 0 . Poor absolute value correlation exists 
at M = 1 . 4 l, although the qualitative effect of angle of attack is well 
predicted. Improved correlation in this case would undoubtedly result if 
the experimental value for the tail-alone side -force coefficient were 
known and substituted in place of the estimated value in equations (12) 
and (13). 


Body -Wing Combination 


Wing-body interaction does not occur in pure sideslip in the case 
where the wing position coincides with the body horizontal plane of symme- 
try other than for small effects of wing thickness. This is true since 
the wing plane then coincides with a body-alone crossflow streamline and 
no change in crossflow takes place. At combined angles and particularly 
for the case where the wing is in an extreme, or tangent, location on a 
circular or nearly circular body, interference forces develop which can 
have significant effects on the directional coefficients of airplanes. 

Body-wing in pure sideslip.- The following expression from slender- 
body theory accounts for the increase in body side force due to the 
addition of a tangent wing to the body: 
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The factor denotes the ratio of increase in side-force coefficient 

referred to body area reference Sg and is plotted as a function of the 
ratio of wing semispan to body radius , r/s^, in figure 11. The term 20 
is the usual slender-body-theory force coefficient for a circular body. 

It was shown in reference 6 that, in general, interference coeffi- 
cients derived from slender -body theory can be applied with acceptable 
accuracy to nonslender configurations and to configurations at supersonic 
Mach numbers. In most cases the effect of Mach number is automatically 
taken into account by the side-force coefficient of the reference compo- 
nents to which each interference coefficient applies. Equation (l4) does 
not take into account the effect of Mach number since here the reference 
component side-force coefficient involved (20) is that of the body alone 
which depends only on sideslip angle. A method is sought by which Mach 
number effects can be taken into account by a modification to equation (l4) 
(strictly valid only for slender bodies at M = l). 

Because of the influence of the forebody on crossflow, the two wing 
panels effectively encounter opposing angles of attack, as indicated in 
the following sketch of a high wing -body configuration. 



These effective angles of attack are proportional to sideslip angle k 0 
and cause a lift on the windward panel and a download on the leeward panel. 
These loads create interference pressures on the body. Wing influence on 
the body is assumed to be confined between Mach helices originating on 
the leading and trailing edges of the junctures of the wing root and body, 
as indicated by the shaded area in sketch (d). 



Sketch (d) 
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It is assumed that the differential pressures acting on the body are 
proportional to those between the wing panels and are given approximately 
in coefficient form by two-dimensional theory as: 


Ap _ 2Ax -2Ax _ 4rcp 

^ Vh* 2 - 1 


(15) 


where k is an undetermined proportionality factor and is assumed a 
constant. This pressure coefficient multiplied by the vertical area pro- 
jection of the interference (shaded) area results in a value proportional 
to interference side force. The value of the projected area is given in 
terms of Mach number, body diameter, and wing root chord length by 
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Multiplying equation (15) hy equation (l6a) cr (l6b) results in a measure 
of the Interference force: 
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It is noted that for c r > nd, the ini erference force (eq. (17b)) 

becomes independent of Mach number and reduces to dependence only on body 
diameter. Division of equation (17b) by (l7e) results in the ratio, rj^, 
a measure of the effect of Mach number on the interference force. This 
ratio times the value of K-g(y\ given by sic nder-body theory (valid at 
Moo = l) then results in the value of K^(w) £ -t Moo > 1. 
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and 

T]^ = 1 for c r > 1 nd (18b) 

The value of r)g is plotted as a function of d)/c r in figure 12. 

The effect of Mach number on the interference coefficient given 

by slender-body theory is thus determined. The side -force coefficient 
increment resulting from addition of a wing to a body becomes 

°%r c *B ‘ | ^ < 19 > 

In order to determine the yawing moment contributed by the interfer- 
ence force, its center of pressure must be known. It is assumed that the 
interference pressures are uniform over the shaded area shown in 
sketch (d). Thus, the centroid of area is identical to the center of 
pressure of the interference force. The longitudinal distance to this 
point from the body nose is 
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The braced terns in these equations give the value of c*/c r or the 
distance to the center of pressure from the juncture of the wing leading 
edge and root chord in terns of root -chor d lengt h, c r . These values are 
presented in figure 13 as a function of (JwJ ' -1 d)/c r . The center-of- 
pressure position being given, the yawing -moment coefficient increment 
can be computed: 
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Body-wing in combined sideslip and angle of attack .- It was pointed 
out that two distinct effects of angle of attack on the directional coef- 
ficients of body -tail configurations occur: the cross -coupling of upwash 

and sidewash velocities, and the nonlinear effect of body vortices. Higher 
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order slender-body solutions which account for body-wing cross -coupling 
are unavailable at present* However, an approximate method is developed 
with which the effects of angle of attack on body-wing interference coef- 
ficients can be taken into account at supersonic Mach numbers. The effect 
of body vortices on body -wing interference is assumed unimportant. 

Consideration will first be given to the midbody wing location. As 
stated previously, the interference coefficient for such combinations is 
negligible in pure sideslip. An important effect arises, however, when 
large combined angles of attack and sideslip are imposed. This occurs 
since viscous crossflow, occurring at large angles of attack, is inhibited 
along a body length equal to the wing root chord. Thus the viscous or 
nonlinear forces acting on the body are reduced. Rough estimates of the 
amount by which the addition of a wing reduces the coefficients contributed 
by the body are 
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The value for ^b(W) depends on wing location and for mid-wing configu- 
rations can be approximated as 7^+(l/2)c r . 

Consideration will now be given to the effect of large angle 
inclination on tangent wing locations. The hypotheses are applied to a 
body with a high wing; however, their application to a body with a low 
wing is obvious. Because of the effect of the forebody on crossflow 
previously discussed, the wing panels are at different angles of attack. 

In addition, the section of the body separating the wing panels is exposed 
to an environment of increased dynamic pressure and reduced Mach number 
prevailing on the compression surface of the lifting wing. These effects 
are represented in sketch (e). The change in effective angle of attack 
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of the panels due to crossflow. Ax, Is again proportional to sideslip, 
or is It is now assumed that the static-pressure field, p, generated 

"by each wing panel in the vicinity of the body is given by two-dimensional 
shock-expansion theory for a flat plate at angle of attack a. This pres- 
sure varies nonlinearly with angle of attack. A pressure difference exists 
between the wing panels since they are at different effective angles of 
attack as indicated in sketch (f). 



Angle of attack 

Sketch (f) 


This pressure differential may be determined from the derivative of p/p^ 
with a along with supersonic flow relationships between M^, q^, and p^ 
as given in reference 15- The differential pressure in coefficient form 
is separated into a linear plus a nonlinear (in square brackets) term 
times the differential angle 2k3. 
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It is to be noted that equation (2k) differs from equation (15) by this 
additional nonlinear term. The resultant reaction force on the body is 
proportional to this differential wing pressure t imes t he combi ned effects 
of local dynamic pressure and Mach number, or q^Hx) 2 -1 / q^vM 2 -1 = rj^, act- 
ing on the area of the interference zone as given by equation (l6) with 
M replacing The following expression can then be written for the 

side -force interference coefficient 
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When a = 0°, M = V^ } q = q^, and = JS 

da /TTp 


to equation ( 17 a) in coefficient form 
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equation ( 25 ) reduce 
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Dividing equation (25) by equation (26) gives a factor k(a) which take, 
m o account the effect of angle of attack on wing -body interference. 
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The ratio k(a) is plotted against angle of attack for various values of 
d/cj. and m figure 14. Taking into account the foregoing effects o^ 
angle of attack and Mach number, the following general expressions result 
for the incremental coefficients. 
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Values for Z B (w) are determined with the aic of figure 15 which gives 
the center-of -pres sure location of the interference pressures of tangent 
wing -body combinations in terms of deviation from the value c, , _ for 
the pure sideslip case given in figure 13 . It is assumed that both the 

pressure and viscous force components of interference act at the same 
location. 
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good moment correlation. This point will be discussed later in conjunction 
with a possible method of increasing the directional stability contribution 
of wing -body interference. 

The experimental and theoretical values for the incremental side- 
force coefficient which result from the addition of a wing in a midbody 
location is presented in figure 16(a). Estimated values are given by 
equation (22) on the premise that the sole effect of the wing is to 
inhibit the nonlinear effects of viscous body crossflow. This simplified 
method somewhat underestimates the incremental coefficient in the higher 
range of angles of attack. The discrepancy indicates a differential load- 
ing between wing panels which is not taken into account. As pointed out 
in the previous section, a differential wing loading acts on the body when 
the wing is asymmetrically located (high or low) on the body. To some 
extent a differential wing loading can also result between sweptback wing 
panels because of unequal sweepback of their leading edges when in an 
attitude of yaw. This effect, which is not taken into account in the 
method, in conjunction with the approximate account which is taken of 
viscous crossflow suppression might explain the lack of better agreement 
between estimate and experiment. 

Experimental and tneoretical values for incremental side-force 
coefficients contributed by wings of various plan forms in high and low 
location on a body arc presented in figures 16(b) through (i). Agreement 
between estimated and experimental values is generally good. This good 
correlation is somewhat surprising since the factor k(cx) in equation (28) 
applies strictly to two-dimensional wings and, in particular, to wings 
with unswept leading edges, whereas the experimental data apply to wing- 
body configurations, among which the wing sweepbacn varies greatly. The 
effect of wing sweepback is twofold: (l) Sweepback reduces the effect 

of k(a) primarily by reducing the effective value of T]^ in equation (27). 
(An approximate estimate by linearized v/ing theory for a sweptback sonic 
leading- edge case indicates a decrease in this term by about 35 percent 
below that of the two-dimensional value.) (2) The unequal panel leading 
occurring with yawed, sweptback wings generates additional differential 
lead on the body, tending to compensate for the reduced value of the fac- 
tor t^. No attempt is made herein to estimate the magnitude of the wing 
sweepback effects. In view of the correlations of figure 16, it appears 
that in lieu of refined methods, equation ( 2'8 ) can be applied within 
±20-perccnt accuracy to combinations having sweptback, but supersonic 
leading edges. 

Correlation of experimental and estimated yawing moments (fig. 17) 
are somewhat inconclusive, since, for the conventional wing locations 
considered here, the center of the wing -body interference pressures is 
close to trie moment reference point, resulting in small values for yawing 
moment. These small yawing moments are possibly no larger than the yaw- 
ing moments of the wing alone (neglected In the method). Fur the mo re, 
as stated in the body-alone discussion, the center of pressure of the body 
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crossflow loading is not accurately known. Thus, exact knowledge of wing 
effects on the body crossflow load distribution is precluded. Estimates 
of the wing-body interference moments by the methods presented is war- 
ranted in cases where the mid-point of the wing root chord is substantially 
displaced from the moment reference point. Previously it was shown that 
large gains in side-force coefficients over that of the body alone occur 
with angle of attack if the wing is located in a high body position. This 
increased side force produces no significant stability increases in the 
cases considered since the force is confined to the vicinity of the wing 
root chord and near the moment center reference. Modifying a wing plan 
form so that the root chord is displaced rearward without changing the 
center-of -pressure location is a means by which wing-body interference 
can be exploited for increased directional stability. Such a wing com- 
pared with a conventional wing in combination with a body is shown: 



interference forces 


Sketch (g) 


Body-Wing-Tail Combinations 


Uncoupled configurations . - The previous sections considered separately 
body-tail configurations and body-wing configurations. For a certain 
class of body -wing -tail configurations, it is only necessary to add the 
wing -body interference terms to the body-tai L coefficients to determine 
the coefficients of a body in combination wi zh a tail and a wing at low 
angles of inclination. The class of aircraft to which this procedure 
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applies is termed the uncoupled wing -tail configuration. An illustration 
of this class is shown in sketch (h). If a body section separates the 



wing from the empennage as indicated so that the Mach lines emanating 
from the wing (containing the wing pressure fields) do not impinge on the 
tail surfaces , the airplane is classed as uncoupled. Actually, an addi- 
tional influence of the wing exists because of changes in sidewash caused 
by the wing addition. The altered sidewash persists to some extent aft 
of the trailing-edge Mach lines. This wing influence on the tail was not 
considered In the method for the sake of simplicity since the accuracy of 
the estimates was not seriously compromised as will be shown. The assump- 
tion of the lack of wing-tail interference applies to the condition of low 
angle of attack and sideslip. Potential flow uncoupling does not preclude 
wing-tail interference in general. When the wing is upstream of the tail, 
the effect of wing vortices on the tail can be significant and should be 
considered. To this end, wing vortex strengths and positions can be 
determined by an obvious extension of the procedure of reference 13 . 
Theoretical body vortex positions, as affected by the presence of the 
wing, are not at present amenable to ready solution for the general case 
of combined sideslip and angle of attack. However, when sideslip angle 
is small, the approximate effect of the wing on body vortex strength and 
displacement can be taken into account by the methods of reference 12. 

Uncoupled configurations in pure sideslip .- The incremental side- 
force and yawing -moment coefficients due to the addition of an upper tail 
to a wing -body configuration are 
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The incremental coefficients of body -wi lg -tail combinations including 
a ventral fin are 
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Uncoupled configurations in combined sideslip and angle of attack . - 
The effect of angle of attack on the incremental coefficient contributions 
of an upper tail to a wing -body combination is taken into account by addi- 
tional cross -coupling and vortex terms. The side-force and yawing -moment 
coefficient contributions of an upper tail to a body -wing combination are 
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The incremental coefficients applying to body -wing -tail configurations 
which include a ventral fin are 
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Closely coupled configurations.- On many aircraft the wing is in 
such proximity to the empennage that the tail surfaces are totally envel- 
oped by the wing flow field. Configurations of this type are classed as 
closely coupled wing -tails, as typified in tie following sketch: 
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Sketch (i) 


Closely coupled configurations in pure sideslip . - In accordance with 
figure 2. the side-force coefficient contribution of an upper tail to a 
body -wing configuration when total coupling exists, is 

< “' y bfv~ Cy bw = k v(bw)C Yv (38) 


The corresponding incremental yawing -moment coefficient is 
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The incremental coefficients pertaining to combinations including a 
ventral fin are 

Cy bwvu‘ c y B w = k v(bw) c y v +k U (bwv) c y u m 

and 
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Values for the K factors are given in figure 18 as a function of body 
radius to tail semispan ratio. This figure gives the interference coef- 
ficients for mid, low, and high body locations for a wing having a ratio 
of semispan to body radius of 5 . More extensive sets of interference 
coefficients are presented in reference 6, wherein the effect of wing 
span and body ellipticity as variables is considered. 

Closely coupled configurations in combined sideslip and angle of 
attack . - Two distinct interference effects between wing and tail are 
caused by angle of attack for closely coupled configurations. One of 
these, the effect of the wing on body vortices, has been pointed out in 
connection with uncoupled wing -tail configurations. The second is the 
change in the side-force coefficient of the tail panel when it is in the 
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wing flow field. This latter effect can he approximately taken into 
account by effectiveness factor discussed in- the body-wing section. 

For a closely coupled wing and tail, the effects of wing vortices are 
precluded since the wing vortex wake is downstream of the empennage. 

The side-force and yawing -moment coefficient increments due to the 
addition of a tail to a body-wing combination in combined angle of attack 
and sideslip are 
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The increments for a body -wing -tail configuration including a ventral 
fin are 
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The effectiveness factor lfy is plotted as a function of angle of attack 
for various Mach numbers in figure 19. The use of this figure is made 
clear by the example on figure 19(a). A suitably scaled outline of an 
airplane with its root chord length coinciding with length c r on the 
chart locates the tail in the proper influence zones. For the example 
shown values of = O .96 for the upper tail and ify = 1.04 for the 
ventral fin are determined by integration of with tail area in each 

influence zone . 
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Partially coupled configurations . - Frequently, body -wing -tail 
configurations have partial coupling between wing and tail. In fact, a 
configuration uncoupled at low Mach numbers may become partially coupled 
at high Mach numbers. Sketch (j) illustrates an example of partial 



coupling. In this class of configurations, the extent of tail area 
affected by the wing depends not only on wing location but on angle of 
attack as well. The method of calculating interference coefficients due 
to the addition of the tail requires only a simple modification to the 
interference coefficient used in the equation which applies to closely 
coupled configurations. The interference coefficient Kv(BW) rno ^“ 

fied as in the method of reference 6. An effective K f is determined 
by the tail area, Sy(BW), lying in the wing expansion (or compression) 
field subject to influence Ky(gy) and that of the tail area, Sv(b)> 
lying out of this field and subject to influence Ky(B) 

K » = K V(BVO S V(BW)+Ky(B)Sv(B) (46) 

s v 

Wing vortex effects are generally negligible for partially coupled 
wing -tail configurations. 

Comparison of estimate and experiment . - Estimated and experimental 
side-force and yawing -moment coefficients of body -wing -tail configurations 
are compared on the basis of the difference between the coefficient value 
for the body -wing -tail combination and that of the body -wing combination. 
This difference is the measure of the tail effectiveness in the presence 
of a winged body. Discussion of the comparison of experiment and estimate 
is confined to side-force coefficient correlation since the degree of 
corresponding yawing -moment correlation follows in close parallel. Esti- 
mated and experimental values are presented in figures 20, 21, and 22. 

Except in one case (fig. 20(d)) the configurations all have a 
forebody extending forward of the wing. Body vorticity therefore affects 
the tail effectiveness and is taken into account in the theory. The 
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effect of the wing on body vortices is taken into account by an approxi- 
mation to the method of reference 1. It is assumed here that both vortex 
strengths and positions are those about a body foreshortened by the length 
of the wing root chord. This approximation then neglects the effects of 
the wing flow field on vortex position; however, the error introduced at 
small sideslip angles is small. Confirmation in this regard is indicated 
in that the correlation of experiment and estimate (figs. 20(a) through (c)) 
is generally good at 3 = 5.1° and poor at p = 10°. Poor correlation 
at the larger sideslip angles can be expected since 10° of sideslip is 
well into the angle range where viscous lateral crossflow effects become 
important. The wing does not suppress and may even augment body vortices 
generated by "strong” sidewash. For instance, data for the low wing con- 
figurations at p = 10° indicate a large difference between experiment 
and estimate. Experimental data for the high wing configuration at 
3 = 10° is, however, in better accord with estimate. The strongest influ- 
ence of the wing on body vortices at large angles of sideslip is thus 
associated with low wing position. 

The wing wake function r^j, based on twc -dimensional shock-expansion 
theory, cannot be expected to apply to wing plan forms of low aspect ratio. 
This consideration was discussed briefly in the section on body-wing inter- 
ference in connection with highly swept wings . A case in point is involved 
in the correlation of experiment and estimate in figure 20(d). It is to 
be noted that the wing plan form is of low aspect ratio and has subsonic 
leading edges at the test Mach number. Application of the charts of fig- 
ure 15 , which apply strictly to two -dimens ioral wings, or approximately 
to wings with supersonic edges, is not valid for this case. A rough cor- 
rection for the effect of plan form was made as follows: The ratio of 

the pressure coefficient along the wing center line (ignoring the presence 
of the body) to that of a two-dimensional wirg at the same angle of attack 
and Mach number was determined by linear theory. The term (l-r^) is 
assumed proportional to the wing pressure coefficient. The magnitude of 
this term is then reduced by the ratio of firite to two-dimensional wing 
center-line pressures yielding an effective \alue for (l-rj^) from which 
the wing wake function T]y is extracted. This correction applies only 

to the wing surface center line since flow at out the finite wing is coni- 
cal. It is therefore of significance that tie corrected estimate gives 
the best correlation for the short span tail whose area is in closer 
"average” proximity to the wing surface than that of the long span tail. 
Another effect of wing plan form is indicated: in the correlation presented 
in figure 20(e) indicating the tail contribution variation with angle of 
attack. The wing plan form in this case was swept forward. Because cf 
the resulting large forward displacement of the wing tips with respect to 
the tail location, tip effects communicate dcwnstream to the tail plane 
and relieve the "two-dimensional” nature of the wing wake field. The full 
compressive wing wake pressures are net then sensed by the (ventral) tail. 
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In the usual case of complete or partial wing-tail coupling, there 
are two effects of wing location on tail side-force contribution: (l) At 

pure sideslip, a wing mounted low on a body diverts crossflow over the 
top of the body which augments the sidewash at the upper tail location. 

This interference effect, taken into account by the interference coeffi- 
cient Ky ( bw)j results in greater incremental tail side force in the 
presence of low wing as compared with high wing cases. This effect of 
wing location is particularly evident in the correlation of estimate and 
experiment presented in figures 21(a) and (b). In all cases, at low angles 
of attack, the tail contributions are substantially greater for the low- 
wing configuration. (2) For a given longitudinal wing location, the wing 
expansion field envelopes a greater portion of the upper tail surface when 
the wing is mounted low on the body. As a consequence, the effectiveness 
of the tail then decreases more rapidly with angle of attack than if it 
were mounted on a high wing -body combination. Thus, at high angles, the 
tail contribution in the presence of a low wing may be no larger than in 
the presence of a high wing-body. This is particularly true at high Mach 
numbers (M^ > 2) as can be verified from figures 2l(a) and (b). Correla- 
tion of estimates with experimental tail contributions at lower Mach num- 
bers are presented in figures 22(a) through (d). At ^ = 1.4-1 the 
configuration of figures 22(c) and (d) falls marginally in the "uncoupled" 
case by the criteria of the method; thus, little effect of wing location 
at a = 0° is predicted. However, experimentally a favorable effect of 
low wing location is clearly evident. This discrepancy between estimate 
and experiment indicates that wing-body effects on sidewash persist to 
some extent downstream of the wing trailing-edge Mach line for marginally 
uncoupled cases. This effect is assumed to wash out downstream of the 
wing trailing edge in the simplified concept presented. However, in view 
of the generally fair correlation between estimate and experiment, the 
simple criterion of coupling and the application of available K factors 
as presented herein and in reference 6 are justified. 

The effects of wing location, of course, apply in the opposite sense 
in regard to ventral tails. The effect of high wing location coupled with 
"minimum" body vortex effects is particularly conducive to increasing 
ventral tail coefficients with angle of attack as shown in figures 20(e) 
and (f). 

It may be concluded from a consideration of the combined effects of 
wing flow fields and body vortices that no optimum tail plan form can be 
specified for maximum directional stability. A short span tail is strongly 
affected by destabilizing body vortices but is adaptable to locations 
least affected by the unfavorable effects of the wing expansion field. 

A large span tail, conversely, is less sensitive to body vorticity but 
usually projects into the unfavorable region of the wing expansion field. 

A compromise must be made on the basis of minimizing the combined effect 
of both destabilizing phenomena at design Mach number. 
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Increased directional stability, provided sideslip angle is not 
large, can be gained with low-aspect -ratio wings whose root chord extends 
to the body nose. Such a body -wing combination generates a wake in which 
the vorticity is concentrated at the wing tips and is therefore "diverted" 
from the directional stabilizer as shown in sketch (k). Stability benefits 

Body vortices (short 

chord wing) 





extended wing 


result even in cases of vanishingly small aspect ratio wings (strakes) as 
noted in references 5 and 16 because of control of body vorticity and 
favorable pressure changes on the forebod; r . 


Horizontal Stabilizers 


Aft located pitch stabilizers. - A derailed consideration of the 
effects of horizontal stabilizers on directional stability of an airplane 
was omitted in the present report. A systematic experimental study of 
horizontal stabilizer effects was not undertaken. A small aft stabilizer , 
simulating an all-movable pitch control, was used on two of the models 
and contributed no important increments in directional stability. Gener- 
ally, this is not the case. Interference effects of a horizontal stabi- 
lizer incorporated in the tail on the directional coefficients of an 
airplane at low angles of inclination can be estimated by the methods of 
reference 6. The interference effects at large angles of inclination 
then can be accounted for by the methods discussed herein in connection 
with closely coupled wing -tail configurat:.ons. In application of the 
method for determining horizontal stabilizer interference effects, certain 
additional considerations must be borne in mind. These considerations 
may require special computations. For instance, the pitch control might 
be operating in the wing flow field and the superposed effects of wing 
and control expansion fields on the vertical tail will be involved. In 
this case the effective control angle can be modified by the wing flow 
field and should be taken into account. 
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The vertical location of an aft mounted pitch stabilizer is affected 
hy optimum longitudinal stability considerations which might not be con- 
sistent with optimum directional stability gains. These two factors must 
be compromised, and specific recommendations in regard to maximizing 
directional stability contributions of an aft located horizontal stabilizer 
are not justified. 

Canard pitch controls.- The interference effects of canard controls 
on the body are identical to those of a wing on the body. The stabilizing 
moment arm for a canard extends forward of the airplane center of gravity 
and thus a negative interference force, as can be attained with a mid-wing 
mounted surface at high angles of attack, contributes positive directional 
stability. Canard locations other than mid-wing contribute positive inter- 
ference forces to the body which cause negative directional stability 
increments . 

The addition of a canard to an airplane modifies the body wake 
affecting the vertical stabilizer. This modified wake can improve the 
stability contribution of twin stabilizers (ref. 17 ) but generally 
decreases the effectiveness of a single stabilizer. 


CONCLUSIONS 


Correlation between estimated and experimental values for side -force 
and yawing -moment coefficients of supersonic aircraft configurations 
presented herein support the following conclusions: 

1. The estimation procedure is sufficiently reliable and consistent 
to enable reasonably accurate prediction to be made of the static stability 
of proposed aircraft configurations. 

2. Lack of good correlation between the present method and experiment 
in certain instances at large sideslip angle is traceable to the effects 

of the wing on body vortices. Further experimental and theoretical 
investigation on body-wing vortices at large sideslip angles is needed 
for a more complete estimation procedure. 

3. The study indicates that significant reduction in the decay of 
directional stability coefficients of supersonic airplanes, generally 
associated with increasing Mach number and angle of attack, can be 
achieved through the use of wings and stabilizers of suitable plan form 
and body location. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., Sept. l6, 1958 
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APPENDIX 

EXAMPLE CALCULATION OF BODY VORTEX EFFECT ON TAIL 


Body vortex effect on the tail will he calculated for the condition 
p = 5°, a = 12°. The center-of-pressure position of the tail at which 
point vortex influence is to he determined is at a distance (x-x s ) /a = 17.6 
from the point of vortex separation. The hody radius ho tail semispan 
ratio is a/sy or r/sy = 0,2. 

The effective angle of attack of the hody is a 1 = Ja 2 + p 2 = 13°. 

At (x-x s ) /a = 17.6 and a 1 = I 3 0 , the correlation presented in reference 11 
gives the approximate strengths and positions for the hody vortices in 
the coordinate syste m Y 1 ,Z * where Z T is in the direction of resultant 
crossflow velocity, s/v^+w 2 . 


L B 


2jta \f 


7 1 Y 1 

= 1.2, — = 1.90, — = O .75 


v^+w^ 


a 


The coordinate positions are then resol\ed into vertical tail 
reference coordinates with cp = tan" 1 (Y ! /Z T ) hy the following transforma- 
tions as indicated in sketch (z). 



Sketch (z) 



H3 *~b 
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s v 

£1 

s v 

— 

s v 

£2 

s v 



The above expressions give for the positions of Pg, hx/sy = 0.007 and 
1 / Sy = 0.413., and h 2 /sy = 0.285 and fa/ s y = 0.294. Corresponding values 
or iy are, respectively, - 1.9 and - 1.2 as determined from figure 8 . 
herefore Ai t = -1.9 - (~1.2) = -0.7. The necessary values are then in 
hand to determine the destabilizing side-force coefficient on the tail 
(>fy(gg) la terms of the side-force coefficient of the isolated tail if at 

an inclination of P = a' = 13 ° as given by equation ( 9 ). 
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TABLE I.- MODEL BODY COOFDINATES 


Station 

B 1 

b iL 

I 


B 

3 

r 

r 

a 

b 

a 

b 

0 

.500 

1.000 

1.500 
2.000 

2.500 
3.000 

0 

0 

.151 

.271 

.351 

.406 

.445 

.474 

0 

0 

0 

0 

3.500 

.151 

.495 

.185 

.123 

.185 

.123 

4.000 

.292 

.508 

.358 

.238 

.358 

.238 

4.500 

.412 

.516 

.505 

.336 

.505 

.336 

5.000 

.481 

.523 

.588 

.392 

.588 

.392 

5.500 

.518 

.529 

.635 

.423 

.635 

.423 

6.000 

• 536 

.536 

.656 

.437 

. 656 

.437 

6.500 

.542 

.542 

.664 

.442 

.664 

. 442 

7.000 

.548 

.548 

.671 

.447 

.630 

.477 

7.500 

• 554 

• 554 

.679 

• 453 

• 574 

• 535 

8.000 

.560 

.560 

.686 

.458 

.536 

. 586 

8.500 

. 566 

.566 

.693 

.462 

.507 

.631 

9.000 

• 570 

• 570 

.698 

. 465 

.486 

.668 

9.500 

• 572 

.572 

.701 

.467 

.472 

.693 

10.000 

• 573 

• 573 

.701 

.468 

.468 

.701 

10.500 

• 572 

.572 

.700 

.468 

.468 

.700 

11.000 

• 570 

.570 

.694 

.468 

.468 

.694 

11.500 

. 566 

.566 

.685 

.468 

.468 

.685 

12.000 

.561 

.561 

.672 

.468 

.468 

.672 

12.500 

• 554 

.544 

.655 

.468 

.468 

.655 

13 . 000 

.545 

.545 

.635 

.468 

.468 

• 639 

13 . 500 

• 535 

• 535 

.611 

.468 

.468 

.611 

14.000 

.522 

.522 

.584 

. 468 

.468 

.584 

14.500 

.508 

.508 

.552 

.468 

.468 

.552 

15.000 

.492 

.492 

.517 

.468 

.468 

.517 

15.625 

.468 

.468 

.468 

.468 

.468 

. 468 
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(a) Uncoupled wing -tail configurations. 
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T 
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BWVU 
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(b) Closely coupled wing -tail configurations. 

Figure 2.- Symbolic definition and order of application of interference 
coefficients K used in determining side force. 













Angle of attack, a, deg 

Figure 5.- Variation of side-force coefficient Cy^ of circular bodies 

with angle of attack. 
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Figure 7*~ Cross -coupling factor Kq> as a function of ratio of body 
radius to tail semispan, r/sy or r/sy. 








AS 



Figure 9*~ ConclulecL. 
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Tail plane 
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"Rectangular" body B A 
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Figure 10.- Vortex formation about a circular and a rectangular body at 
combined angles of attack and sideslip (a = 8°, P = 5°). 

















Figure 13.- Center-of-pressure location, c, of wing -induced side force 
on tody in terms of wing-root chords fron juncture of wing leading 
edge and root chord, a = 0°. 







Angle of attack, a, deg 


(a) High wing location. 

Figure 14.- Interference coefficient k(a) as a function of angle of 
attack for various Mach numbers . 











Angle of attack, a, deg 

(t>) Low wing location. 


Figure 15 . - Concluded 





















(a) Mid-wing location. 

Figure 18 .- Interference coefficients K for tails in presence of an 

r/s^ =0.2 body -wing combination. 



Figure 18.- Continusd. 
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(c) High wing location 
Figure 18 . - Concluded 


(b) a - 10’, M„ = 2 

Figure 19.- Charts of tail effectiveness parameter 







irjBBP wnnr.fl 



1 

J 

.48 


T 

t 

mm 


mttKmnwm 

R9CSH 



SSSSafl 



E lgiBBBBBBBBBglBBBBBal 

«Ea .-. ■ ..■■ •.•■ .= ,'.: ■ - ■" • ■-'■ -■■•. -81 

a BaMWBS«l : x m . im 

MmSS 

























Angle of attack, a, deg 
(c) 


Angle of attack, a, deg 
(d) 


Figure 20.- Comparison of experimental and estimated side-force and 
yawing -moment coefficient increments resulting from addition of 
vertical tail to various body -wing combinations at = 2.9^. 







Angle of attack, a, deg Angle of attack, a, deg 

(a) High wing (b) Low wing 


Figure 21.- Comparison of experimental and estimated side-force and 
yawing -moment coefficients resulting from addition of vertical tail 
to a body at various Mach numbers; unpublished data from tests in the 
Langley Unitary Plan wind tunnel. 
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Figure 22.- Comparison of experimental and estimated side-force and 
yawing-moment coefficient increments resulting from addition of 
vertical tails to tody -wing combinations at ^=1.41 and 2.01. 
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